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Abstract 

A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed 
by solving numerically the coupled differential equations of motion for inflaton and photon mode 
functions from the onset of inflation to the end of reheating. The backreaction due to particle 
production is also included self-consistently. We find that this particular inflation model realizes 
the idea of a warm inflation in which a steady thermal bath is established by the particle production. 
In most cases this thermal bath exceeds the amount of radiation released in the reheating process. 
In the strong coupling regime, the transition from the inflationary to the radiation-dominated 
phase does not involve either a preheating or reheating process. In addition, energy density peaks 
produced near the end of inflation may lead to the formation of primordial black holes. 
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I. INTRODUCTION 


The inflation model, in which our observable Universe is only a tiny local patch of a 
causally connected region that underwent an exponential or de Sitter expansion at early 
times driven by an inflaton potential, is generally accepted for explaining the observed 
spatially flat and homogeneous Universe. A simple version of the model such as the slow- 
roll inflation predicts quasi de Sitter vacuum fluctuations during inflation which could give 
rise to Gaussian and nearly scale-free metric perturbation containing both matter density 
fluctuations (scalar mode) and gravitational waves (tensor mode) [I|. 

Recent astrophysical and cosmological observations such as cosmic microwave background 
anisotropies, gravitational lensing, baryon acoustic oscillations, supernovae, and Hubble 
constant have supported the slow-roll scenario, and further constrained the non-Gaussianity 
of the matter perturbation and the tensor-to-scalar ratio [^. This enables us to discriminate 
between different types of inflation or even disfavor certain models j^. However, some 
fundamental questions have yet to be answered. What is the nature of the inflaton? Do 
classical matter density inhomogeneities in the present Universe genuinely come from vacuum 
fluctuations of the inflaton? Are the matter density fluctuations truly Gaussian? Is the 
tensor-to-scalar ratio predictable? Future cosmic microwave background measurements and 
large-scale-structure surveys should be able to address some of these questions or pose a 
challenge to the standard slow-roll model. 

There has been a lot of studies on inflationary models that go beyond the simplest 
single-field, slow-roll inflation. Recently, a class of models has considered a new source for 
generating metric perturbation during the standard slow-roll inflation through a coupling 
between the inflaton and other quantum fields, with particular attentions to the influences 
of particle production and its associated backreaction to the slow-roll inflation. This leads 
to very interesting effects such as the suppression of large-scale density fluctuations j^, 
the generation of non-Gaussian and non-scale-invariant density power spectrum the 

damping of primordial gravitational waves due to a copious production of free-streaming 
particles [8|, the generation of stochastic gravitational waves that could be detected at 


loll, and the formation of primordial black 
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ground-based gravity-wave interferometers [Oj, 
holes at density peaks near the end of inflation 

In particular, a lot of attentions have been paid to the study of the phenomenological 
effects on inflation via a pseudoscalar-vector coupling 0,0 M3, 


a 




( 1 ) 


where / is an energy scale and a is a dimensionless coupling constant. The pseudoscalar (p is 
an axion-like held that drives inflation, — is the vector held strength tensor, 

and = ^e^'^°‘^Fai3 is its dual. The pseudoscalar can be a pseudo Nambu-Goldstone boson 
in the natural inflation 171 or the many of its variant models [l^, which is expected to couple 
to some gauge held as the one shown in Eq. ([T]). In addition, axion inflation is generically 
at work in a broad class of inflationary models in supergravity, in which reheating after the 
end of inflation requires considering the pseudoscalar-vector coupling [l^. So far, tackling 
the inflaton-vector system has relied on an adiabatic approximation under which simple 
analytic solutions for the vector mode equations can be used. Although this approximation 
is valid in the beginning of inflation, it becomes unreliable near the end of inflation and 
even inapplicable in the era of reheating, when physical effects of main concern are taking 
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place. Recently, the authors in Ref. [19| have studied using lattice simulations the onset of 


the reheating epoch at the end of inflation. They found that for a wide range of parameters 
preheating is efficient and that in certain cases the inflaton transfers all its energy to the 
vector fields within a few inflaton oscillations. In this paper, we perform a numerical analysis 
of the system from the onset of inflation towards the end of reheating, solving self-consistently 
the coupled equations of motion of the fields including field perturbation and backreaction. 
In particular, we pay attention to the transition from the end of inflaton to the preheating. 
We will review the the inflation model in Sec. Hi] and set up the numerical scheme in Sec. IIIIl 
Sec. IIVI is our numerical results that are compared to those in Ref. 
conclusions in Sec. El 


19 . We make our 


II. PSEUDOSCALAR INFLATION 

We consider a simple model of inflation driven by a pseudoscalar that couples to a U{1) 
gauge field via the interaction ([1]). The action is given by 


s= 


1 1 a ~ 

-fR- - U((p) - - -(p F^^ F,. 


( 2 ) 


where R is the curvature scalar, Mp is the reduced Planck mass, V{(p) is the inflaton poten¬ 
tial, — dyAp is the field strength tensor, and F^^ = Fa ^/is its dual. 

Note that l/y/—g is added to the dual tensor because is a tensor density of weight —1. 
Here we assume a spatially flat Friedmann-Robertson-Walker metric: 


ds'^ = —Qpydx^dx^ = a^{rj){drf — dx^), 


(3) 


where a{ri) is the cosmic scale factor and rj is the conformal time related to the cosmic time 
by dt = a{r])dT]. The Hubble parameter H = d/a has a conformal time analogue R = a'/a, 
where the dot and the prime denote derivatives with respect to t and g, respectively. 

From the action ([2]), we can write down the Friedmann equation, the equation of motion 
for the inflaton, and the Maxwell equations, respectively: 


= 


3M2 




+ 2R^ _ vV + a' = a^-E • F, 

orj^ or] dp j 


d^A 


a dp ^ 


a. 


dA 


- + V(V ■ A) = -ttV X a - -Vp X —, 


drf 

|-(V . 1) = . (V X i"), 

or] f 


f dr] 


f 


dr] 


(4) 

(5) 

( 6 ) 
(7) 


where for the Maxwell equations we have chosen the temporal gauge, i.e. Ap = (0, A), and 
we have introduced the physical “electric” and “magnetic” fields. 


F = - 


1 dA 
a? dr] ’ 


B = —V X A. 
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( 8 ) 













The above equations form a complete set of coupled differential equations. In the present 
consideration, we study the production of gauge quanta by the rolling inflaton via the 
interaction during a slow-roll inflation, taking into account self-consistently the backreaction 
of the gauge quanta production on inflation. Furthermore, we consider the generation of 
curvature perturbation in this scenario. 


A. Gauge quanta production and backreaction 


To calculate the production of gauge quanta, we separate the inflaton into a mean held 
and its huctuations: 

(p = 0(t) + (5(p(t, f). (9) 

To the hrst order in quantum helds, we can consistently impose the condition, V ■ A = 0. 
Then, Eqs. (E]) and ([7]) become a mo diked wave equation. 


dr]'^ 


- V^A- X A 


0 . 


To proceed, we decompose the gauge held A{ri,x) as 


Av,x) = 

\=± 


d^k 

(27r)3/2 


ex{k)axik)Ax{T], + h.c. 


( 10 ) 


( 11 ) 


where the annihilation and creation operators obey 


aA(fc),ap?) = SxyS(k-k'), 


( 12 ) 


ex are normalized circular polarization vectors satisfying k ■ e± yvj = 0, k x e± yvj = 

^ike± (^kj, e± (— {k^ , and ex (pj ■ (p = Inserting the decomposition flTT|) 
into Eq. flTOl) . we obtain the equation of motion for the mode functions. 


'P. 

dr]'^ 


+ k‘^ T 2aHk^ 


A±{r],k) 


0 , 


a<p 

Wh 


(13) 


It is well known that either one of the two modes exhibits a spinoidal instability as long as 
the modes satisfy the condition, k/{aH) < 2|,^|. When the inhaton rolls down the potential, 
these unstable modes grow exponentially by consuming the inhaton kinetic energy. As such, 
the energy density of the produced gauge quanta is given by the vacuum expectation value. 


h{E^ + B^) = ^ldkk^'£A'x\'‘ + k^\Ax\^). ( 14 ) 

A=i 

Also, we have 
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On the other hand, the production of gauge quanta gives rise to a backreaction on the 
background. The background evolution is therefore governed by the mean held parts of 
Eqs. dl]) and ([5]) with the backreaction effects included, i..e 




3H^ = 


M2 




( 16 ) 

( 17 ) 


B. Inflaton fluctuations and curvature perturbation 


To study the generation of inflaton fluctuations during inflation, one can integrate out 
the gauge fields in the real-time formalism to obtain an effective equation of motion for 
the inflaton fluctuations, including backreaction and dissipation effects. Similar work about 
integrating out scalar fields has been done in Refs. [3, 1^. However, in Refs. 0, [Tsl it 
was proposed that one can estimate these effects in the regime of strong backreaction by 
incorporating a dissipation term in the equation of motion: 


+3,3/14- 

^ '1 ^ 


9f2 


dt 


d02 


(X / —y —y —* —* 

= - (E-B-{E-B 


(18) 


where 


/? = 1 - 27r4 


a{E-B) 


(19) 


/ 3H<j) 

Furthermore, it was shown that the solution of this equation can be well approximated 
by 0,0 


Sip = 


a 


(e-b-{e-b 


( 20 ) 


3l3fH^ 

which leads to a contribution to the power spectrum of the curvature perturbation given by 
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= {cixY) = -^(V) = 


a{E-B) 


( 21 ) 


When the backreaction becomes large, the second term in Eq. flT^ dominates and we can 
have approximately 




( 22 ) 


III. NUMERICAL SCHEME 


We employ the chaotic inflation scenario [2^ to unravel the particle production and 


its backreaction by full numerical calculations. In such a universe, the standard slow-roll 
inflationary background is driven by the potential 


V(<p) 



( 23 ) 
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where the parameter m = 1.8 x 10^^ GeV characterizes the mass of the inflaton field. Here we 
rescale all dynamical variables in terms of the reduced Planck mass, Mp = 2.435 x 10^® GeV. 
Hence, m = 7.39 x 10“®. We set / = 1 and 0o = —14.15 and 0o = 6 x 10“® respectively 
for the initial position and the initial velocity of the inflaton (p. The e-folding since the 
beginning of inflation is dehned by N{t) = Jq where it is chosen that Uq = 1. 

The background evolution (f){t) and a{t) with feedback from gauge quanta production can 
be obtained numerically by solving self-consistently Eqs. (H, (USD, (USD, and fll7p , where the 
mode functions of the gauge held during inhation are governed by Eq. flT^ . with the initial 
conditions of the gauge mode functions given by A±q = \l\p2k and A±o = —ik/y/^. To 
solve this set of coupled differential equations, we use the second-order Runge-Kutta method. 
To accurately capture the feedback of the gauge quanta production in Eqs. ffTSP and flTTH . we 
need to integrate over the range of k to which the A± modes have contributions at a given 
time. Initially, all A± fc-modes have small contributions. However, the H+ modes develop 
spinoidal instability and grow exponentially when they satisfy the condition k/{aH) < 2^. 
On the other hand, the backreaction of the fc-modes decay with the growth of a, so one 
should only count modes that satisfy (8^)“^ < k/{aH) < 2^ [^. Note that the range is 
time varying during inflation. However, we include all fc-modes with 10“^ < k. The scale 
factor a will grow up to about at the end of inflation, so we have to calculate a very 
large range of fc-modes (10“^ < k < 10^^). Since it is difficult to calculate very high-fc modes 
that are varying too fast with time, to reduce the computing load we approximate these 
modes by their initial values until they enter the range with k/{aH) < 2^. Although this 
would cause phase shifts of the mode functions, it can hardly affect the time evolution of the 
mode amplitudes. Since only the mode amplitudes are used in the background equations, 
the approximation works very well as far as the background evolution is concerned. Figured] 
shows the range of /c-modes that are included in the backreaction at a given time during 
inflation. 

To calculate the backreaction, we need to integrate all relevant fc-modes by the Trapezoid 
rule, the Simpson’s rule, and/or the Boole’s rule. In practice, because of the range of k 
relevant to the backreaction is too big, it is problematic for us just to divide the whole range 
into equal integration intervals. As such, we do the following trick. Initially, the range of k 
is narrow and the value of k is small. During inflation, the range of k becomes larger and 
larger, so we choose different partitions for different fc-modes. The partition in low-/c region 
is small whereas that in high-fc region is wide. In practice, we adopt equal partitions in 
\og-k range. By integrating the modes in each interval and then summing all intervals we 
thus obtain the total backreaction effect to the background evolution. Furthermore, we will 
increase the time steps as well as the fc-resolution until we obtain steady results. 

Because of limited computational resource, in some cases we have tested our concluding 
results in toy inflation models that have shortened duration of inflation while capturing rel¬ 
evant features of the particle production and the backreaction to the background evolution. 
In this toy inflation model, the range of k is considerably reduced and therefore we can save 
more computational power to increase the time steps, fc-resolution, and fc-range for calcu¬ 
lating the growth of the fc-modes and their feedback. We have found that this strategy is 
very useful to eliminate computational artifacts and ascertain the physical results obtained 
in realistic inflation models. 
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FIG. 1: Solid line is the upper limit and the dashed line is the lower limit. The shaded area 
between the two limits represents the range of fe-modes at a given e-folding N that contribute to 
the backreaction. k in this figure and all dynamical variables in the following figures are rescaled 
by the reduced Planck mass Mp = 2.435 x 10^® GeV. 


IV. NUMERICAL RESULTS 


Let us first consider a strong coupling a = 32. In this case, ^ = 2.2 in the beginning of 
inflation at V = 0 and inflation ends aX N = 64. This value of ^ is the upper bound coming 
from measurements of the power spectrum and the bispectrum of the cosmic microwave 
background j^. The influence of the backreaction of the gauge quanta production to the 
inflaton motion is shown in Fig. ([2]). Because the range of fc-modes for the growing gauge 
held is related to 0 that in turn affects the backreaction, the temporal variation of 0 is 
very complicated during inhation. This is the reason that the scalar held 0 and the Hubble 
parameter H exhibit bouncing at the late stage of inhation. We have also conhrmed the 
occurrence of this bouncing behaviour in a toy inhation model with the strategy as mentioned 
in Sec. IIIIl The numerical results for the evolution of 0 and H are shown in Fig. [3] and 
Fig. m respectively. It is apparent that the backreaction slows down the inhaton motion 
and thus prolong the duration of inhation by about 10 e-foldings. It is also interesting to 
see that inhation ends gracefully without preheating or reheating process. These results are 
consistent with those found in Ref. [l^. Here we show detailed features of the evolution 
towards the end of inhation. Figure |5] shows the evolution of the parameter The fast 
oscillatory behaviour of ^ at the late stage of inhation puts a caution that one cannot use 
the analytic solution for the growing mode H_|_, which has assumed an adiabatic condition 


that treats ^ as a constant parameter (see Ref. [12| and therein previous work), for the whole 
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course of inflation, especially near the end of inflation. 

The results of the power spectra of the curvature perturbation in Eq. (|2T]) for various 
a parameters are plotted in Fig. [6l where each power spectrum has already included the 
amplitude of the vacuum inflationary perturbation as predicted in a standard slow-roll in¬ 
flation, which is given by Note that we have adjusted the value of the initial 0o 

in each case such that the duration of inflation is 64 e-foldings. In general the results in the 
hgure are consistent with the large coupling approximation for fl22D . The effect of the 
backreaction can modify drastically the power spectra found in Ref. [l2| near the end of in¬ 
flation. The dashed line in Fig. [6] is the upper bound on the power spectrum coming from the 
non-detection of primordial black holes (referring to Fig. 5 of Ref. 12| and noting that the 


bound ends at about N = 54). We hnd that the black hole bound constrains a < 20, which 
is consistent with, though slightly tighter than, the Endings in Refs. 12|, ll3|. For a > 20, 
the power spectra show prominent peaks that exceed the black hole bound. Interestingly, 
we observe that near the end of inflation all power spectra show high spikes that seem to be 
beyond the black hole bound when presumably extrapolated to A > 54. Note that in some 
of the curves, the power spectrum spikes above unity and is cut off there. We warn that in 
the present consideration the power spectrum is no longer valid when its value becomes near 
unity. In this situation, one must take into account the effects of gravitation, which could 
presumably dampen the spikes to a consistent level j2l|. However, the occurrence of these 
density spikes suggests that primordial black holes may probably be formed at high-density 
regions near the end of inflation. If so, a more stringent bound on the value of a can be 
derived. However, more studies about the formation of black holes at the density spikes in 
the present consideration are needed to confirm this suggestion. 

In Fig. [21 the dashed line is the result without backreaction. When the term {a/f){E■ B) 
on the right-hand side of Eq. flTB]) is negligible, 0 ~ 0. The slow-roll condition that 0 
keeps steady until the end of inflation is satisfied. Because the initial conditions are 0o < 0 
and 00 > 0, the backreaction {a/f){E ■ B) <0. The backreaction kicks in when either 
one of A±{k) grows exponentially for k satisfying the growth condition k < acpa/f. Once 
the backreaction starts to affect the motion of the inflaton, the inflaton decelerates, i.e. 
0 < 0. Consequently, 0 is reduced and therefore the largest growing mode of A±{k) in the 
backreaction stops growing. Then, the backreaction {a/f){E ■ B) of the produced gauge 
quanta will be redshifted as a“^, thus allowing 0 to build up itself and letting more k- 
modes grow exponentially in the backreaction to affect the inflaton motion again. Since the 
backreaction effect of the gauge modes is not immediate, even if a mode starts to affect the 
inflaton motion, the k value of the mode could be much smaller than the growth threshold, 
a(j)a/f. It will then keep exponentially growing until 0 is reduced below the growth condition. 
As such, there occurs a delay on the response of the backreaction to 0. This explains why 
the power spectra show spikes in Fig. O 

Now we consider the reheating temperature of the pseudoscalar inflation. The cou¬ 
pling ([1]) implies that an inflaton would decay into two gauge bosons. The perturbative 
decorate, i.e. when the mean held 0 = 0 or the background is the ground state, is given 
by 0 

9 

_ am 

^ 64^' ^ ^ 

When the inflaton condensate decays, the Universe is reheated with the gauge quanta. The 
gauge quanta may not be thermalized; however, we treat them as a thermal bath with a 
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FIG. 2: Phase diagram of the inflaton mean field (p and its velocity (p, with the initial (pQ = —14.15, 
cpo = Qx 10“®, and the coupling constant a = 32. This case corresponds to ^ = 2.2 in the beginning 
of inflation {N = 0). The solid line does take the backreaction into account, while the dashed line 
does not. 


temperature defined by: 

+ = (25) 

where Qes is the number of effectively relativistic degrees of freedom at the time of reheating. 
Here we assume Qes = 106.75. This perturbative reheating temperature can be estimated 
by comparing the decay rate with the Hubble rate, T ~ 3H: 


Tper ~ — 


Stt 


/_m V''" (K\ 
W) [mJ \ f ) 


Mf ~ 4.42 X 10“'”a. 


(26) 


However, the spinoidal instability of the gauge held in the presence of the background inhaton 
mean held may lead to a striking reheating process. Figure [7] shows the temperature of the 
gauge quanta in Eq. fl2^ from the beginning of inhation to the time of reheating. Apparently, 
a higher value of the coupling constant a would lead to a higher reheating temperature. In 
the case of a = 32, because of strong backreaction, in general cp is smaller than that in 
the standard slow-roll inhation, thus resulting in a longer duration of inhation. There is 
not a preheating phase; namely, (p gradually reduces to zero, rather than undergoing a 
period of held oscillations. For small coupling constants with a < 1.45, the backreaction 
is no longer important because of the small growth of gauge mode functions for all fc’s. 
In these cases, to evaluate the fc-integration in Eq. fflTl) . we can approximate the mode 
functions by their initial amplitudes, x4± 1/y/^ and A'^ ~ /c/2, and integrate over the 

range 0 < /c < 2^aH. Hence, we obtain + B‘^)/2 ~ 2^‘^H‘^/tt'^. Then, the reheating 
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FIG. 3: Evolution of (/> as a function of the e-folding N in the case with a = 32. The solid line 
does take the backreaction into account, while the dashed line does not. 



FIG. 4: Time evolution of the Hubble parameter H in the case with a = 32. The solid line does 
take the backreaction into account, while the dashed line does not. 
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FIG. 5; Time evolution of the parameter ^ in Eq. (|13p in the case with a = 32. 


temperature can be estimated by 

Tre ~ O.Ua^/f ~ 0.84 x lO’^a. (27) 

This is consistent with the numerical calculations of the reheating temperatures for a < 1.45 
in Fig. [3 Note that the production of the gauge quanta during inflation can sustain a 
nearly steady temperature flTT)) which is much higher than the perturbative reheating tem¬ 
perature fl26l) . The pseudoscalar inflation with an axion-photon-like coupling is a concrete 
example that realizes the idea of the existence of a thermal bath during inflation in the 
so-called warm inflation jl^. In most warm inflation models an ad hoc frictional term (T^) 
reflecting the decay of the inflaton is added to the inflaton equation of motion and the re¬ 
heat!^ process is through the perturbative decay of inflatons. Also see, for examples in 
Ref. [2^, some effort to derive the frictional term from hrst principles. 

In the present work, we have solved the coupled differential equations numerically with 
backreaction included self-consistently. Our numerical method is quite similar to that 
adopted in a previous paper that considered the same axion inflation model but in a differ¬ 
ent context of magnetogenesis [2l[. In the paper, the authors pointed out that for strong 
couplings the generation of inflaton fluctuations S^p by the gauge quanta would become too 
large to invalidate the linear approximation, and therefore they put forth a consistency con¬ 
dition for linearity, where = 0^/2 -|- R(0) is the inflaton energy density. 

Hence, they found that a < 8.4 and that for larger a the non-linear coupling between 
and the gauge quanta would require performing full lattice simulations. 

However, in Ref. [2l[, the inflaton perturbation equation used to calculate the inflaton 
fluctuations sourced by the gauge quanta does not include a dissipation term. This dis- 
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FIG. 6: Time evolution of the power spectra in Eq. (I2ip during inflation. The lines shown in the 
figure, in order from left to right, are for a = 32,28,23,20,16, respectively. These correspond in 
turn to ^(N = 0) = 2.2,1.91,1.53,1.32,1.04. Note that we have adjusted the value of the initial c^o 
in each case such that the duration of inflation is 64 e-foldings. The black hole bound is the upper 
dashed line extended to about N = 54. 


sipation term becomes significant in the strong coupling regime, as we have discussed in 
Sec. IIIBI When a becomes large, the inflaton fluctuations d(p in Eq. ffTK]l are damped by 
dissipation and their variance is given by = p/(jia)^ as shown in Eq. fl22D . Moreover, 

the inflaton fluctuations are sourced by the gauge quanta that come from the inflaton energy. 
Therefore, we propose to improve the consistency condition to that rn?‘{5p) < ptot, where 
Ptot = is the total energy density. Let us denote the energy density of gauge quanta 

by pem = {E'^ + B‘^)/2. Then, we have ptot = + Pem- During inflation, ptot — P<p and the 
condition is reduced to w?{5ip) < p^p. During reheating, 6ip can be generated by the gauge 
quanta even though pp pem and the condition becomes rn?{5p) < pem- 


Now we compare our numerical results with the lattice simulations made in Ref. [19 


Figure [H] shows the time evolution of the ratio pEiu/Ptot from the onset of inflation to the 
end of reheating for a wide range of values for a. Note that inflation ends at p = 0 when 
t = tend- Also, before inflation the time is reckoned by e-foldings AN = N — A^(tend) and 
after inflation it is conformal time dehned by p = /a{t') with rescaled a(tend) = 1- 

In Ref. 19|, the authors simulated the preheating and reheating processes by tracing the 
evolution of pEn/ptot with time after the end of inflation (p > 0) for a ranging from 7 to 
13. We compare our results for the ratios to those in Ref. [l^ with the common input 
parameters and the same time ranges. We hnd that overall the two results agree with each 
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FIG. 7: Time evolution of the gauge quanta temperature T defined by /2 = (tt^ /30)5effT^ 

where = 106.75. The lines shown in the hgure, in order from top to bottom, are for a = 
32,14.5,1.45,0.145,1.45 x 10“^, 1.45 x 10“^, 1.45 x 10““^, respectively. The left black dot on each 
line indicates the end of inflation timed by —H/H^ = 1. The right black dot denotes the beginning 
of the radiation dominated era, after which T oc a~^. 


other quite well while there are differences in detailed sub-structures. Here we are mainly 
concerned with the inflaton dynamics as well as the gauge particle production near the end 
of inflation. It is worthy to further scrutinize the preheating and reheating processes and 
compare the two methods in more details. 

To check the consistency condition for linearity, we have calculated using Eq. fl2ip 

and plotted the time evolution of the ratio w?{5ip^)/ptot in Fig. [9l The results in the hgure 
are consistent with the large coupling approximation for fl2^ . It seems that in some 

of the curves the ratio exceeds unity and violates the consistency condition. In fact, before 
the ratio reaches unity, the reheating process has already hnished as shown in Fig. |8l For 
instance, the ratio for the curve with a = 9 reaches unity at the time r/ ~ 7 while Fig. [8] 
shows that the reheating for this case is hnished at r/ ~ 5. 


V. CONCLUSIONS 

We have constructed a numerical code for computing self-consistently the motion of 
inhaton, the particle production, and the density perturbation in the pseudoscalar inhation 
model with the axion-photon-like coupling, aipFF/{4:Mp), where Mp is the reduced Planck 
mass. We have used a quadratic inhation as a working example and investigated the ehects 
of the coupling for a < 32. This model is unique in the sense that the photon held receives 
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AN rjim ^) 

FIG. 8: Time evolution of the ratio of the energy density of gauge quanta, /?em = 
and the total energy density, ptot = Inflation ends at /S.N = 0 or t/ = 0. Before inflation 

the time is reckoned by how many efoldings prior to the end of inflation. After inflation the scale 
factor is rescaled to unity and the conformal time p is in units of m~^. The lines shown in the 
figure, in order from bottom to top, are for a = 7,8, 9,11,14,16, 23,32, respectively. 


an effective mass that is proportional to the inflaton velocity. In particnlar, that one helicity 
state of the photon obtains a negative mass leads to a spinoidal instability. As a consequence, 
a wide spectrum of photons is produced during inflation efficiently enough to form a thermal 
bath that leads to a reheating temperature, Tre ~ which is much higher than the 

perturbative reheating temperature, Tper ~ Furthermore, we have calculated 

the density power spectrum induced by the backreaction of the photon production on the 
inflaton fluctuations. Using the bound on the density power spectrum coming from the 
formation of primordial black holes, we have derived a more stringent upper limit a < 20. 
In addition, a strong coupling can produce high density spikes in the density power spectrum 
near the end of inflation that may form primordial black holes. 
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FIG. 9: Time evolution of the ratio of the inflaton fluctuations, rn?{5ip^), and the total energy 
density, ptot = The lines shown in the figure, in order from bottom to top, are for 

a = 7,9,14,16, 23, 32, respectively. 
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